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Abstract. We present the first detection of CO emission lines in the Ho- filaments at distances as far as 50 kpc from the centre 
of the galaxy NGC 1275. This gas is probably dense (> 10 3 cm~ 3 ). However, it is not possible to accurately determine the density 
and the kinetic temperature of this relatively warm gas (T kin ~ 20-500 K) with the current data only. The amount of molecular 
gas in the filaments is large — 10 9 M (assuming a Galactic N(H2)/Ico ratio). This is 10% of the total mass of molecular gas 
detected in this cD galaxy. This gas has large-scale velocities comparable to those seen in Ha. The origin of the filaments is still 
unclear, but their formation is very likely linked to the AGN positive feedback (Revaz et al., 2008) that regulates the cooling 
of the surrounding X-ray-emitting gas as suggested by numerical simulations. We also present high-resolution spectra of the 
galaxy core. The spatial characteristics of the double-peaked profile suggest that the molecular web of filaments and streamers 
penetrates down to radii of less than 2 kpc from the central AGN and eventually feed the galaxy nucleus. The mass of gas inside 
the very central region is ~ 10 9 M Q , and is similar to the mass of molecular gas found in the filaments. 

Key words, galaxies: cooling flows — galaxies: individual: NGC 1275 — galaxies: ISM — galaxies: kinematics and dynamics 
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1. Introduction 

The cD galaxy NGC 1275 lies at the centre of the Perseus clus- 
ter, the brightest galaxy cluster in the sky in X-rays. The galaxy 
is at aredshiftof 0.01756, (D A = 72.6 Mpc; 1" is 350 pc). Past 
detections of CO in NGC 1275 were mainly toward the centre 
of the galaxy (Lazareff et al., 1989; Mirabel et al, 1989; Reuter 
et al., 1993; Braine et al., 1995; Inoue et al, 1996; Bridges & 
Irwin, 1998; Lim et al., 2008). We also studied CO in filamen- 
tary structures all around the galaxy (inside r<10 kpc) and in 
three pointings at galactocentric distances as large as r~25 kpc 
(Salome et al., 2006, 2008a, 2008b). The molecular gas coin- 
cides with the Her filaments (e.g., Cowie et al. 1983; Conselice 
et al. 2001), within the hot gas seen in soft X-rays at 0.5 keV 
(Fabian et al. 2006). These results were in good agreement with 
the CO observations of the central 15 kpc made by Reuter et 
al., (1993) although the regions observed were not exactly the 
same. The central region was also observed with the VLA by 
Jaffe (1990) who detected extended 21cm absorption in front 
of the radio source 3C84. The HI is elongated in a east-west 
direction like the CO and the inferred neutral hydrogen mass 
found is ~5 10 y M o , that is, of the same order as the molecular 
gas mass (~10 10 M o ) found in the same region from CO obser- 
vations (Salome et al. 2006). Irwin & Bridges (2001) found a 
north-west/south-east extended continuum emission at 450yt/m 
and 850/mi with SCUBA. This core-subtracted emission excess 



reveals a large amount of cold dust inside the central ~20 kpc 
(6xl0 7 M Q with T dust ~20K and/3=1.3). 

In this paper, we describe a new search with the IRAM 30m 
telescope for molecular gas in the Ha filaments at the largest 
radii from the centre of NGC 1275 (outside the already ob- 
served central region and as far as 50 kpc). We clearly identify 
molecular filaments in five regions (see Fig. 1): 

( 1) A 30 kpc-long filament, running north-south, not previously 
observed in CO (regions rl through r4). The largest radius at 
which we detect CO (region rl) is 50 kpc (in projection) north 
of the centre of the galaxy. 

(2) The southeast Ha filaments (rlO, rl 1, rl5). The distance 
between the end of this filamentary system (region rl 1) and the 
end of the north filament (rl) is 200" (70 kpc), 

(3) A tangential filament to the north-east (rl9, r20), where we 
have extended our previous CO observations. 

(4) The northwest area, including the well-known horseshoe 
filament, was observed in regions r5, r6, r8, r9, and rl6. 

(5) The central region around 3C84 was re-observed. Here, we 
made a high-sensitivity spectrum with high spectral resolution. 

Section 2 presents the observations. Section 3 describes 
the new CO detections. Section 4 discusses the excitation con- 
ditions of the molecular gas at large radii around NGC 1275 
and the possible origin of the very distant molecular filaments. 
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2. Observations 

We used the IRAM 30 m telescope on Pico Veleta, near 
Granada, Spain. Observations were made in December 2007. 
Receivers were tuned to the CO(l-O) and CO(2-l) lines, red- 
shifted to the adopted systemic velocity of NGC 1275 (1 13.280 
and 226.559 GHz). We used two 512x1 MHz filters at 3 mm 
and two 250 x 4 MHz filters for the 1.3 mm receivers, thereby 
covering 1300km/s at both wavelengths. We used wobbler 
switching, and checked the pointing on 3C84, the radio source 
at the centre of NGC 1275. We observed in excellent weather 
with T sys typically 250 K at 1 .3 mm. To study the large extent of 
molecular gas around NGC 1275, we observed 18 regions in- 
side the Ha filaments (Fig[TJ- The beams of the 30 m telescope, 
at 3 mm and 1.3mm, are 22" and 11" respectively (correspond- 
ing to 7.7 and 3.8 kpc at NGC 1275). The line intensities are 
given in main-beam brightness temperatures, T m b- The ratio of 
T m b to T* at the 30 m is T m b/T* = F e g/B e ff = 1.27 at 3 mm and 
1.75 at 1.3 mm (see http://www.iram.es/). 

The data were calibrated with the MIRA software and re- 
duced with the CLASS package. We dropped spiky channels 
and bad scans, and subtracted linear baselines from each spec- 
trum before averaging. Table [T]summarises the results. 



3. Results, by region 

3.1. The gas masses 

For comparison with other work, the gas mass was estimated 
from a standard Milky Way conversion factor M gas /Lp =4.6 
M (K.km/s.pc 2 ) -1 (Solomon et al. 1997) for both transitions. 
Note however that the ICM metallicity is known to be less 
than solar Z<0.6Z o in the Perseus cluster core (Schmidt et al., 
2002). This low metallicity would lead to a larger N(H2)/CO 
conversion factor than the standard one used here (Leroy 
et al, 2009) and thus under-estimated gas masses. The CO 
luminosity is defined by L' CQ = Tb AV Q.s D^ with Tb Qs the 
source brightness temperature times the source solid angle, AV 
the line width and the angular distance of the source. This 
gives L' CQ = 23.5 Q s »b Ico Dj^ (1+z), with z the source redshift, 
Q s »b the solid angle of the source convolved with the telescope 
beam and Ico=J T m b AV the CO intensity (T m b being the main 
beam temperature). This is for L' CQ in K.km/s.pc 2 , Q, s *b in 
arcsec 2 , Da in Mpc, and Ico in K.km/s. With this conversion 
factor, we find molecular gas masses of ~3-16xl0 7 M in each 
pointed region of the filaments. We used the same conversion 
factor for both transitions since line ratios were close to one. 
We find more mass in the (larger) 3mm beam than in the 1mm 
beam. So the CO emission may be extended. Note however 
that this is just a hint because of the large uncertainties on the 
line width and main beam temperature measurements for both 
transitions. 



3.2. The 30-kpc-long north filament, at very large 
radius: 

This filament extends over a radial distance from 20 to 50 kpc 
from the galaxy's centre. We detected CO(2-l) at all positions 
observed, but no CO(1-0) apart from r4 (Table [1] regions rl 
to r4). The fact that CO(2-l) is easier to detect means that the 
CO(1-0) is strongly affected by beam dilution, and therefore 
the emission is fragmented and clumpy on a scale of a few 
kpc (in region r4, where both lines are detected, the CO(2- 
l)/CO(l-l) main beam temperature ratio is ~2). The CO(2-l) 
line centre changes from -184km/s (at rl), to -65 km/s (at r2), 
then jumps to +100 km/s (r3), and +50 km/s (r4). These line 
velocities show the same trend as in Ha (Hatch et al. 2006). 
Figure |4] (middle panel) shows the observed position- velocity 
diagram along this filament. There is a positive gradient from 
-200 km/s to +100 km/s as the radius decreases from 50 to 
35 kpc, after which the velocity falls back to 40 km/s at a ra- 
dius of 25 kpc. The good agreement with the Ha suggests that 
the CO traces most of the mass in the filaments, possibly sur- 
rounded by ionised gas radiating in Ha, at the interface with 
the hot intra-cluster medium. The shape of the gradient and 
the velocity reversal at the end could be due to the projec- 
tion of the curved filament. Judging from the other loops that 
are oriented more face-on, it is obvious that many of the fila- 
ments have a significant curvature. The linewidth of regions rl 
to r3 is 200 km/s while region r4 shows a narrower linewidth 
of 100 km/s, more typical of the sub-structures in the eastern 
filament (see Salome et al. 2008b). With the standard Milky 
Way conversion factor described above, we find a molecular 
gas mass of ~3-14xl0 7 M in the northern filament. 



3.3. The bright southeast filaments: 

Regions rlO, rl 1 and rl5 (the 'blue-loop') cover the southeast 
Ha filaments. Regions rlO and rl 1 are only detected in CO(2- 
1), while region rl5 is also detected in CO(1-0). This filament 
has negative velocities only. There is a velocity change from 
-130km/s to -250 km/s between regions rl 1 and rlO (compat- 
ible with the velocities of a falling filament). Region rl 1 has a 
broad CO(2-l) linewidth (370km/s), while in region rlO, the 
CO(2-l) profile is narrow (120 km/s). Slightly farther east is 
the region rl5 at a velocity of - 130 km/s. These southeast Ha 
filaments are wider than the Northern ones, and the CO signal 
is probably less affected by beam dilution. 



3.4. The north-east tangential filament: 

This filament is not radially extended like most of the other 
filaments around NGC 1275, but is instead perpendicular to a 
line through the centre of the galaxy. Region rl9 is detected in 
CO(1-0) and region r20 in both CO(1-0) and CO(2-l). The 
line velocities are all negative, with no clear gradient. The line 
widths are narrower (70, 100 and 180 km/s) than those of the 
other detections. 
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Fig. 1. CO spectra superposed on the Hff filament structure around NGC 1275 (Conselice et al., 2001). The double circles show 
the 1.3 and 3 mm beamwidths at the regions where CO was detected. CO(l-O) and CO(2-l) spectra are shown when detected. 
The channel width is 42 km/s and the y-axis is main beam brightness temperature, in mK. The velocity scale is from -560 to 
560 km/s in each insert. Small circles (1.3 mm beamwidth) show regions where no CO was detected. See Table[T]for more details. 
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Table 1. CO detections and tentative detections, or 3<x upper limits. See Salome et al. (2008a) for complementary data. 



Site Offsets 

"X" 


CO 
line 


T,„h 
mK 


V 
km/s 


AV 
km/s 


^co 
Kkm/s 


1O 8 M 


Northern Filament 


rl -21.0x134.9 
rl 


CO(2-l) 
CO(1-0) 


2.5 ± 1.2 
< 3x0.6 


-184.7 h 


:28.2 


198.9 = 


i 52.9 


0.5 ± 


0.1 


0.3 


r2 -19.7x112.9 
r2 


CO(2-l) 
CO(1-0) 


3.5 ± 1.2 

< 3x1 


-65.3 ± 


23.4 


246.4 = 


= 43.4 


0.9 ± 


0.2 


0.5 


r3 -18.6x90.9 
r3 


CO(2-l) 
CO(1-0) 


2.6 ± 1.2 
< 3x1 


100.2 : 


t24 


187.6 


±48 


0.5 ± 


0.1 


0.3 


r4 -18.6x68.9 
r4 


CO(2-l) 
CO(1-0) 


5.7 ± 1.5 
2.7 ± 0.8 


51.3 ± 
55.2 + 


12.6 
19.4 


112.9 = 
209.3 = 


= 23.3 
= 38.8 


0.7 ± 
0.6 ± 


0.1 
0.1 


0.4 
1.4 


South East Filament 


rlO 28.8X-34.9 
rlO 


CO(2-l) 
CO(1-0) 


13.6 ±2 

< 3x2 


-255.8 


±8 


114.5 = 


= 23.2 


1.7 ± 


0.2 


1 


rll 39.8X-54.2 
rll 


CO(2-l) 
CO(1-0) 


6.9+ 1.7 
< 3x1.8 


-137.5 h 


: 19.4 


373.1 = 


= 47.6 


2.7 ± 


0.3 


1.6 


rl5 55.2X-40.8 
rl5 


CO(2-l) 
CO(1-0) 


7.5 ± 2.4 

6.6 ± 1.7 


-250.6 h 
-241.4 d 


: 21.6 
: 15.5 


146.8 = 
123.5 = 


= 67.2 
= 30.7 


1.2 ± 
0.9 ± 


0.4 
0.2 


0.7 

2 


North East Tangential Filament 


rl9 35.9x15.5 
rl9 


CO(2-l) 
CO(1-0) 


< 3x1.7 
8.9 ± 1.3 


-77.7 d 


:5.4 


68.8 ± 


10.7 


0.7 =fc 


0.1 


1.5 


r20 16.7x25.7 
r20 


CO(2-l) 
CO(1-0) 


7.7 ± 1.7 
6.4 ± 1.6 


-79.8 ± 
-100.4 h 


14.6 

: 12.2 


184.7 = 
109.4 = 


= 28.6 
= 24.3 


1.5 ± 
0.7 ± 


0.2 
0.2 


0.9 
1.7 


Horseshoe Filament 


r5 -28.1x27.3 
r5 


CO(2-l) 
CO(1-0) 


9.1 ±2.3 
< 3x1.6 


23.5 ± 


13.6 


101.5 = 


= 58.0 


1.0 ± 


0.3 


0.6 


r6 -45.6x41.4 
r6 


CO(2-l) 
CO(1-0) 


< 3x2.7 
< 3x2 
















r8 -75.3 26.5 
r8 


CO(2-l) 
CO(1-0) 


3.9 ± 1.5 
< 3x0.9 


-17.1 ± 


18.9 


122.1 = 


= 35.2 


0.5 ± 


0.1 


0.3 


r9 -80.8x40.7 
r9 


CO(2-l) 
CO(1-0) 


< 3x1.6 

< 3x0.9 
















rl6 -52.5x55.2 
rl6 


CO(2-l) 
CO(1-0) 


< 3x1.5 

< 3x1.2 
















Southern Filament 


rl2 0.8X-60.1 
rl2 


CO(2-l) 
CO(1-0) 


< 3x1.7 

< 3x1.6 
















rl3 -15.9X-30.9 
rl3 


CO(2-l) 
CO(1-0) 


< 3x1.8 

< 3x1.3 
















r21 -0.6X-21.1 
r21 


CO(2-l) 
CO(1-0) 


4.7 ± 2.3 
< 3x3.2 


86.7 ± 


37.2 


222 ± 


78.3 


1.1 ± 


0.4 


0.6 


Western Filament 


rl7 -58.0X-10.1 
rl7 


CO(2-l) 
CO(1-0) 


< 3x1.8 
< 3x1 
















Center 


rO 0x0 


CO(2-l) 


72 ± 6.9 


4± 


4 


293 = 


= 10 


22.5 =1 


:0.6 


15.3 



T„,;, = main-beam temperature at line peak; V = velocity relative to systemic; AV = linewidth (FWHM); / C o = integrated line intensity; 
M g as(H+He) is for M/L' co = 4.6 M Q (K kms~' pc 2 )~' (Solomon et al., 1997). Note that the integration time is not the same for all regions, so 
the rms (computed for a spectral resolution of 42 km/s) varies from one pointing to another. 
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Fig. 3. CO(2-l) at the position of 3C84. The axes are main- 
beam brightness temperature (in mK) vs. velocity, relative to 
226.559 GHz (z = 0.01756). The channel width is 5 kms 1 . 



3.5. The horseshoe: the faintest region: 

The horseshoe was re-observed in different regions (r5, r6, r8, 
r9, rl6) all along both sides of the U-shaped loop. We clearly 
detect the CO(2-l) line in the trunk of the filament, but this is 
the only firm detection. CO(2-l) is tentatively detected in re- 
gion r8. The other regions were not detected in CO(1-0) or 
CO(2-l). Surprisingly, this region has weaker CO emission 
than the other filaments. If the inclination of the horseshoe fila- 
ment along the line of sight is large, then the regions we observe 
could be at galactocentric distances as far as the faint regions 
detected in the Northern filament (50 kpc). In such a case, the 
small amount of molecular gas would not be visible with the 
rms noise we reached for those regions. The only exception is 
in region r8 where the integration time was longer (see Table 
1). 

3.6. Gas in the central 4 kpc: 

Figure [3] shows a double-peaked CO(2-l) line at the position 
of 3C84 in 5 kms 1 channels in a 4-kpc diameter region. The 
negative velocities peak at -45 kms -1 and the positive ones at 
+65 km s~' This feature can be interpreted as (i) inflowing gas, 
(ii) outflowing gas or (iii) a rotating pattern as observed in lines 
in the NIR by Wilman et al. (2005). It is not possible to dis- 
entangle from this three possibilities with the present obser- 
vations only. However, in previous observations (Salome et al., 
2006), the absence of a larger-scale rotation pattern was pointed 
out (inside a 8-kpc diameter region). 

The double-peaked line cannot therefore be unambiguously 
interpreted as evidence for a central rotating disk. Instead, the 
large accumulation of molecular gas (nearly 1O 9 M ) in the 
central 4-kpc diameter of the galaxy could be the result of the 
recycling of cold gas inflow from the filaments back into the 
central potential well over several Gyr. 

4. Discussion 

4.1. Large velocity dispersions 

Observations by Salome et al. (2008a) showed CO linewidths 
of 10 to 30 kms -1 in the eastern filament when observed with 
a 3" beam. In our 22" beam, however, the CO(1-0) linewidth 



in that same filament is close to 100 km s _1 . In the current data 
from the 30 m telescope, most of the CO(1-0) and (2-1) lines 
have widths of order 200 km s _1 in beams of 4 to 8 kpc. In the 
northern filament, the velocity dispersion is even larger. These 
large linewidths are difficult to explain by motions within a sin- 
gle filament. It is more likely that in our beam, we observe an 
ensemble of very thin sub-filaments (as observed in the optical 
with the HST by Fabian et al. 2008) that consist of clumps of 
molecular clouds. Furthermore, these individual sub-filaments 
are very likely to have different projection angles, which would 
also increase the observed linewidths, when observed with a 
larger beam. Nevertheless, region rl 1 has a particularly broad 
line width (373 km/s). More observations with a larger band- 
with coverage, as provided now by EMIR (Eight Mixer re- 
ceiver^ on the IRAM 30m-telescope, would be useful to con- 
strain with more confidence the effect of baseline fluctuation 
to which the results are very sensitive in the low signal/noise 
regime. 

4.2. Comparison with numerical models 

In Figure [4] we compare the properties of the CO in the north- 
ern filament with the prediction of models in which the molec- 
ular gas forms out of uplifted ambient gas trapped by a plasma 
bubble. We have modelled the AGN feedback by instanta- 
neously generating bubbles of gas in the ICM as already done 
by several authors (see Briiggen & Kaiser 2001, for instance). 
We have used very-high resolution TreeSPH simulations com- 
bined with a multiphase model and a model of plasma bubbles. 
We have also taken into account the presence of cooler and 
denser around the radio bubbles, assuming that it is the result 
of ambient gas pushed out by the inflating bubble. In our sim- 
ulations, a fraction of the 1 -2 ke V gas present at the centre of 
clusters is trapped and entrained by the rising buoyant bubble 
to higher radius where the AGN heating is less efficient. The 
radiative cooling then makes it cool, thereby forming cold fil- 
amentary structures in the wake and in the rim of the bubbles 
(see Revaz et al. 2008 for a detailed description of the model). 
We tried to fit the data presented here by the model number 2 
of Revaz et al. (2008). This is at t = 260 Myr and for a filament 
that is inclined by 60° with respect to the line of sight (the top 
of the filament pointing towards the observer). The model re- 
produces not only the filament elongation, but also the mass 
of molecular gas detected. At radii < 40 kpc, the mean veloc- 
ities and the velocity dispersions are also in good agreement, 
but at larger radius, the model fails to reproduce the negative 
velocities and most of the time it underestimates the velocity 
dispersions by a factor of 1.5, 2 and 1.3 for regions rl, r2, and 
r3 respectively. The negative velocities can be understood if the 
gas is still being uplifted by the rising bubble. In the simulation 
however, at temperature below 10 6 K, the cooling gas decou- 
ples from the ambient gas and falls towards the centre. This 
discrepancy would be removed if another physical process, like 
viscosity or magnetic fields (Fabian et al., 2008) increase the 
coupling between the different gas phases. 



1 http://www.iram.es/IRAMES/mainWiki/EmirforAstronomers 
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Fig. 2. Results of escape probability models from the RADEX program (Van der Tak et al., 2006) for four different kinetic 
temperatures, a range of H2 densities and a range of CO molecular column density per km/s. The grey scale mapping and the 
black contour lines show the expected CO(2-l) emission T«(co2i) [K] (see wedge to the right of each panel). Dashed blue lines 
show the CO(2-1)/CO(1-0) line temperature ratios 0.6, 0.9 and 1 as predicted by the RADEX program. Dash-dot red lines show 
the predicted CO(3-2)/CO(2-l) line temperature ratios 0.6, 0.9 and 1. The measured CO(2-1)/CO(1-0) line temperature ratio 
in rl5 and r20 are close to one. So the possible solutions lie along the dashed-blue line of value one. We see that it is thus 
difficult to constrain the range of possible kinetic temperatures and H2 densities. Nervertheless, even with a high Ty,, of 500 K, 
the density must be >10 25 . Note also that if we consider a standard Milky Way M g . ds /L' co conversion factor, then we expect 
Nco/AV>l 10 17 cm^km.s -1 if f2D ~10~ 3 ie Tr ~10K, which gives slightly higher densities (>10 29 ). Finally, Bridges & Irwin 
(1998) reported a CO(3-2)/CO(2-l) line ratio close to one in the central ~ 8kpc region. If this is also true for the filaments, 
then the possible solutions would lie at the intersection of the dashed blue line (values ~0.9-l) with the dash-dot red line (values 
~0.9-l), which means even higher densities. The high-J CO lines are thus important diagnostics to determine the gas properties. 



4.3. Line ratios 

Among the 12 regions detected here, only three have been de- 
tected in CO(1-0) and CO(2-l). All the others are detected in 
CO(2-l) but not in CO(1-0). One region is detected in CO(l- 
0) and not in CO(2-l). As already mentioned, this could be 
explained by effects of beam dilution that makes the CO(1-0) 
emission from a clumpy molecular gas more difficult to de- 
tect than the CO(2-l) emission in general. Nevertheless, when 
closer to the centre, the filaments are widespread and overlap 
both the 3mm and 1mm beams. However, the emission seen 
in each wavelength may not come from the same part of the 
filaments. Evidence for that is the different line shapes seen 
in region r4, indicating CO(1-0) and CO(2-l) emission from 
dynamically and spatially different regions. This also explains 
the clear CO(1-0) detection alone in region rl9: the region that 
dominates the emission lies in the CO(1-0) beam but not in the 
CO(2-l) one. In conclusion, the CO(2-l) to (1-0) line ratios 
are very uncertain, as are the parameters derived from these 
values. 

The three regions with best signal-to-noise ratios (r4, rl5, and 
r20) were detected in both CO(1-0) and (2-1). Averaged over 
these three regions, the Ico and T m b (2-1/1-0) ratios are both 
equal to 1 .5+0.3. Because interferometer maps (Salome et al 
2008b; Lim et al. 2008) show the CO is extended in the di- 
rection of the filaments, the different beam areas at CO(1-0) 
and (2-1) mean that the apparent line ratios should be corrected 
downward by a factor of 2, so that the corrected Ico and T m b 
ratios are close to 0.8, that is, of the order of unity. 



Simulations with the line escape probability program RADEX 
(van der Tak et al. 2007), indicate that these observed CO(2- 

1) /CO(1-0) line ratios of order unity can be obtained for gas 
kinetic temperatures ranging from 20 to 500 K, for H2 gas den- 
sities >10 2 - 5-3 crrT 3 andNco/AV >10 15 cnT 2 km/s (right panel 
of Fig. 2). So this gives very poor constraints on the gas prop- 
erties. 

With the JCMT, Bridges & Irwin (1998) observed CO(3- 

2) emission from the central 21-arcsec region of NGC 1275 
and marginally detected the I3 CO(2-l) inside a 22" beam. 
They made a two-component LVG model to reproduce 
the 12 CO(3-2)/ 12 CO(2-l)=1.25+0.25 and 12 CO(2-l)/ 12 CO(l- 
0)=0.74+0.11 line ratios. Their best model is for a cold com- 
ponent (10K) at a density of 10 3 crrT 3 and a hot (170K) compo- 
nent at a density of 3.1 10 4 cirT 3 . Furthermore, new interferom- 
eter observations at CO(3-2) (Salome et al, in preparation) in- 
dicate that the CO(3-2)/CO(2-l) ratio in some central filaments 
is also of order unity. This indicates that the CO is certainly 
optically thick in the central region and that the low-J CO line 
brightness and excitation temperatures are all close to the gas 
kinetic temperature. 

It is not possible to definitively rule out an optically thin 
gas in the outer filament without CO(3-2) observations there. 
However, given that the observed (beam diluted) CO(2-l) 
brightness temperatures are about 7 mK, the area filling factor 
f2D in the CO(2- 1 ) beam must be < 1 0~ 3 . As a check on this rea- 
soning, for the three filament regions (r4, rl5 and r20), the H2 
gas masses (column 8 of Table 1) estimated from the observed 



P. Salome et al.: A very extended molecular web around NGC 1275 



GO 



50 



40 



Qj 

5 30 

CO 



20 



10 



1 1 1 


1 1 


rllj 


M = 1.6 


r2 M 


L M = 0.6 




. M = 0.7 


- 1 




r4 ^ 


N M = 0.8 






1 
1 

1 1 1 


► 

1 1 




-10 -5 5 10 
Position [kpc] 



-400 -200 200 400 
Velocity [km/s] 



n i r 

FWHM = 130 



J L 




rl 



i r 

FWHM = 118 



r2 



J L 



u 



~i 1 r 

FWHM = 140 



J L 




r3 



~l I 
FWHM = 117 



J L 



_L 




r4 



-400 -200 200 400 

Velocity [km/s] 



Fig. 4. Northern filament: comparison of the model with the observations. Left: contours are the surface density of molecular gas, 
predicted from the model. CO(2-l) beamwidths and positions are shown by black circles. Gas masses in each beam, predicted 
from the model, are in units of 10 8 M o . Middle: position-velocity diagram. Pink and blue isophotes in the color map are the 
model predictions. Black crosses are CO(2-l) observations. Horizontal bars are the FWHM linewidths. Red diamonds are Ha 
data (Hatch et al. 2006). Right: velocity distribution of the gas extracted from each region of the model. The predicted FWHM 
linewidths of the distributions are in km/s. 



CO line intensities with a standard Milky Way conversion fac- 
tor, are about 10 8 Mo, within the beam. The beam diameter at 
CO(2-l) corresponds to 4kpc, or 1.3 10 22 cm, so the equivalent 
smoothed-out H2 density within a 4-kpc diameter sphere con- 
taining 10 8 MO of H2 gas would be 0.03 crrr 3 , far too low a gas 
density to collisionally excite the CO (the CO(1-0) and CO(2- 
1) critical densities being 3 10 3 cnr 3 and 10 4 cnT 3 respectively, 
for a temperature of 100K). 

Therefore, to reach the gas density of >10 3 crrr 3 needed to 
make the low-J CO lines optically thick, and to give line ratios 
close to unity, the volume filling factor must be <3 10~ 5 . 
The corresponding area filling factor within the beam must be 
<10~ 3 , and the linear filling factor fiD=f2D = ^3D > P er P en dicular 
to the filaments (or along the line of sight), must be <0.03, rel- 
ative to the 4-kpc beam size. 

These upper limits on the filling factors, derived from the CO 
conversion factor and the gas density required for excitation, 



are thus perfectly consistent with the filling factor upper limits 
derived above, from the observed line temperatures and those 
predicted by the escape probability programs, for the observed 
low-J line ratios. The linear filling factor of <0.03 implies a CO 
filament width (or CO clump size) of < 1 20/ V^V pc, N being the 
number of clouds inside the beam, which is also consistent with 
our size upper limit of <450 pc estimated directly from our in- 
terferometric maps (Salome et al. 2008b). 
There is therefore a coherent picture in which to interpret the 
observed low-J CO line ratios. The key is the low filling factor, 
which explains the observed main-beam brightness tempera- 
tures, the high gas density needed for collisional excitation in 
the escape probability predictions, the fact that the low-J line 
ratios are of order unity, the fact that the low-J lines are op- 
tically thick, and the interferometer upper limits on the CO 
clump sizes. 
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4.4. Molecular filament formation 

The origin of the huge molecular filaments is unknown. We 
discuss here two possibilities. Scenario (1): The molecular gas 
forms in situ, far from the galaxy's centre, from uplifted warm 
gas that cooled down behind rising bubbles, and eventually falls 
back (Revaz et al., 2008). This scenario must explain how to 
make the IT, which is difficult to form without dust (Ferland 
et al. 2009). Certainly, the bubble morphology of the filaments 
and X-ray cavities show that this gas must be mixed from the 
centre of the galaxy up to 100 kpc outward, and some dust from 
the galaxy must have enriched the whole region. The presence 
of dust could augment H2 formation rates even in the presence 
of X-rays. But the sputtering due to hot ions in the X-ray emit- 
ting gas can destroys unshielded dust in < 1 Myr (Dwek & 
Arendt, 1992), so if there is dust in the filaments, it must have 
come from the galaxy's inner regions, rather from the intra- 
cluster medium. Scenario (2): The molecular gas comes from 
the inner part of the galaxy, entrained and dragged out by ris- 
ing bubbles. The problem in this scenario is to explain how 
very small and dense clouds can be uplifted inside warmer gas. 
This is difficult to achieve, even in a highly viscous outflow 
and there is no such example in any observed galaxy yet. More 
modelling is needed to define the limits of these different sce- 
narios. Stronger limits on the amount of dust in the filaments 
would help constrain the first scenario, and modelling the mo- 
mentum transfer to clumps in a viscous flow (with and without 
magnetic viscosity) would help constrain the second scenario. 

4.5. Star Formation 

There is no evidence of strong star formation in NGC 1275 fil- 
aments. Canning et al. (2010) have identified, with the Hubble 
Space Telescope, only two stellar regions in the south and 
south-east filaments (22 kpc from the nucleus). The predom- 
inantly young stellar population is concentrated in the Ho- fila- 
ments. The authors derived a lower limit of 9xl0 8 M o stars in 
blue clusters and a maximum lifetime of ~0.75xl0 8 yr, which 
gives a very large SFR of ~20 M /yr. We have compared the 
amount of molecular gas with the star formation activity in 
these filaments The region of interest is indicated as rl5 ('Blue- 
loop' filament) in the present paper and fills the CO(1-0) beam 
size. We have found a molecular gas mass of ~2xlO 8 M . So 
if the SFR is as large as ~20 M /yr then, the amount of gas is 
about 20% of the amount of recently formed stars. This means 
that most of the gas must have already been transformed into 
stars during the filament dynamical time (that roughly corre- 
sponds to the age of the stars). So if the star formation is still 
taking place in rl5, the star formation efficiency inside this fil- 
ament is very high. Indeed, even if the amount of molecular 
gas has been underestimated because of the sub-solar metallic- 
ity in the filaments Z<0.6ZO, (Schmidt et al., 2002), the fila- 
ments will still appear to be locally the place of an enhanced 
star formation activity with a star formation efficiency closer to 
ULIRGs than to normal spiral galaxies. The gas uplifted by the 
AGN bubble expansion in the ICM is pushed and compressed. 
In these overdense regions, the cooling is more efficient and 
molecular gas forms. This reservoir could then be converted 



very efficiently into stars in places where large scale dynami- 
cal effects can enhance the star formation (like in spiral arms 
of normal galaxies). The interface between the AGN cavities 
and the intracluster medium seems to be the place of such pro- 
cesses. 

5. Conclusions 

We have detected CO in a very extended filamentary network 
around NGC 1275. Up to about 10 9 M o of molecular gas is de- 
tected inside thin and elongated filaments at distances between 
7.8 kpc and 50 kpc from the galaxy centre. It is possible that 
the massive filaments seen in CO are an effect of the AGN 
feedback. The large scale motions in the Northern filament are 
difficult to explain with free-fall models. An accurate study of 
the dynamics of the filament should take into account the vis- 
cosity of the warmer gas (some part of it being uplifted). This 
will help lower the model predictions of the cold gas veloci- 
ties. These single dish observations do not resolve the very thin 
filaments as observed with the HST (Fabian et al. 2008), so 
we certainly also mixed several filaments and averaged their 
individual velocities. The CO observed linewidths are large. 
This confirms that the emission very likely comes from sub- 
structures unresolved by the 30 m telescope, possibly molecu- 
lar cloud complexes within the filaments. Finally we reported a 
velocity structure in the emitting gas within 2 kpc of the AGN 
that could be a rotation disk. Whether this is the fate of the 
molecular gas found in very distant regions (20-50 kpc) is still 
an open question. 
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